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The VERITAS lACT observatory has carried out an extensive survey of the Cygnus region between 67 and 82 
degrees in galactic longitude and between -1 and 4 degrees in galactic latitude. This region is a natural choice 
for a Very High Energy (VHE) gamma-ray survey in the Northern Hemisphere, as it contains a substantial 
number of potential VHE gamma-ray emitters such as supernova remnants, pulsar wind nebulae, high-mass 
X-ray binaries, and massive star clusters, in addition to a few previously detected VHE gamma-ray sources. It 
is also home to a number of GeV gamma-ray sources, including no less than four new high-significance sources 
detected in the first six months of Fermi data. The VERITAS survey, comprising more than 140 hours of 
observations, reaches an average VHE point-source flux sensitivity of better than 4% of the Crab Nebula flux at 
energies above 200 GeV. Here we report on preliminary results from this survey, including two source detections, 
and discuss the prospects for further studies that would exploit the joint coverage provided by VERITAS and 
Fermi data in this region. 



1. Introduction 

To date, only a few moderate-scale surveys have 
been performed in gamma rays between 100 GeV and 
lOTeV: the HESS scan of the central region of the 
Galactic Plane the much less sensitive HEGRA sur- 
vey of the quarter of the Galactic Plane between —2° 
and 85° in galactic longitude [l|, and the VERITAS 
scan of the Cygnus region under discussion here. How- 
ever, there is a strong motivation for performing such 
surveys; an unbiased search of a substantial region of 
sky is less subject to the experimental and theoretical 
prejudices that guide most VHE gamma-ray observa- 
tions, and therefore (as demonstrated by HESS) offers 
greater scope for serendipitous discoveries. An unbi- 
ased search also allows for more quantitative state- 
ments to be made about the source population in the 
region surveyed. 

The stereoscopic imaging atmospheric Cherenkov 
telescope (lACT) array VERITAS has just completed 
a two-year survey of a 5 by 15 degree portion of the 
Cygnus region of the Galactic Plane. The Cygnus re- 
gion was a natural target for survey observations, as 
it is already known to contain a significant number 
of potential TeV gamma-ray emitters. In the GeV 
(20MeV - 300 GeV) energy band, it is home to a 
number of sources or potential sources, including no 
less than 4 distinct Fermi sources Moreover, both 
Fermi and its predecessor EGRET have detected dif- 
fuse emission from this region that is greater in flux 
than all of the currently resolved sources taken to- 
gether. Viewed in the energy range between lTeV-50 
TeV, it contains a pair of unidentified TeV sources 
(MGRO J2031+41 and MGRO 2019-^37) detected 
by the Milagro Gamma Ray Observatory, a water 
Cherenkov extensive air shower array that is sensitive 
to TeV sources at a median energy of 20 TeV over 
the entire sky Q, as well as the unidentified source 
TeV J2032-t-4130 (first detected by the HEGRA lACT 
array [l^), that is spatially coincident with MGRO 
J2031+41. The exact nature of these sources is cur- 



rently unknown. There is also a significant catalog 
of objects detected at other wavelengths, including 
SNRs, pulsar wind nebulae (PWNe), high-mass x-ray 
binaries (HMXBs) and massive star clusters, that are 
considered potential TeV sources. 



2. Survey Observations 

Survey observations, which began in April 2007 and 
were completed in November 2008, cover the field be- 
tween galactic longitudes 67° and 82° and galactic 
latitudes — 1 ° and 4° with a grid of pointed observa- 
tions. Grid points have 0.8 ° separation in galactic lon- 
gitude and 1 ° separation in galactic latitude, allowing 
for substantial overlap in the fields of view for observa- 
tions at nearby grid points. Approximately one hour 
of observing time was taken at every grid point (gen- 
erally within a 1-3 day period) , with that hour broken 
into 20-minute observation periods. Figure [T] shows 
both a schematic of the survey observation strategy 
and the acceptance-corrected ("effective") exposure 
time for the "base" survey (that is. the survey data 
taken prior to follow-up observations) over the entire 
survey field. The base survey achieves a relatively uni- 
form effective exposure of 6 — 7 hours. The full sur- 
vey dataset (> 140 hours of good-quality observation 
time) has regions of enhanced exposure due to follow- 
up observations taken in Fall 2008, Spring 2009, and 
Fall 2009. 

Survey data have been quality-selected to remove 
runs with unstable trigger rates, poor weather condi- 
tions, and known hardware problems. Almost all sur- 
vey data were taken on moon-less nights, with a few 
percent taken under slight to moderate moonlight con- 
ditions. Since survey observations began during the 
commissioning period for the VERITAS array, data 
in the base survey were taken in two different config- 
urations. Observations from Spring 2007 were taken 
with a three-telescope array configuration (including 
telescopes 1, 2, and 3) and observations from Fall 2007 
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Figure 1: EfTective exposure map for the VERITAS Cygnus sky survey before follow-up observations, based on data 
obtained from Spring 2007 through Fall 2008. The color scale to the right gives the effective exposure time in minutes. 
The black box indicates the boundary of the survey region proper; the black circles show examples of the overlapping 
fields of observation that are tiled to produce the survey. 



and Spring and Fall 2008 were taken with the full four- 
telescope array. 

Survey observations were taken over a range of 
zenith angles between 10° to 35°. Observations were 
scheduled in such a way as to keep the average zenith 
angle of observations for any survey pointing at 20°, 
but constraints in terms of available time did not al- 
ways allow for this. We estimate that 45% of the 
survey region is covered by observations taken at an 
average zenith angle of 20° or smaller, another 33.5% 
by observations taken at 20° < z < 25°, 9% by ob- 
servations taken at 25° < z < 30°, and 12.5% by 
observations taken at z > 30°. 



2.1. Follow-up Observations 



Follow-up data taken in Fall 2008 were also taken 
with the four-telescope array configuration: observa- 
tions during this period are a combination of two types 
of observations: wobble observations at 0.5° to 0.7° 
offsets from points of interest (> 4(7 fluctuations or lo- 
cations of reported AGILE or INTEGRAL transients) 
within the survey region, and additional exposure on 
survey pointings near a point or region of interest. 
Further follow-up observations ( 60 — 70 hours total) 
using a wobble strategy were taken in Spring and Fall 
2009. Due to a scheduled VERITAS upgrade, data 
taken in Spring 2009 were taken with a three-telescope 
array configuration containing telescopes 2, 3, and 4, 
while data taken in Fall 2009 was taken with the full 
array in an upgraded configuration, where telescope 1 
has been moved to improve the instrument baselines. 
This upgrade has improved the sensitivity of the array 
to weak point sources by about 20% 0. 



3. Survey Analysis 

Given that the H.E.S.S. survey revealed a popula- 
tion of Galactic TeV gamma-ray emitters that was 
biased towards hard-spectrum, moderately extended 
sources 0, and the fact that that the Milagro survey 
of the Cygnus region 0] showed sources of significant 
apparent extent, it was reasonable to expect that some 
or all of the visible gamma-ray sources in this region 
would also have relatively hard spectra and significant 
extension. 

In order to optimize the survey's sensitivity to 
such sources while not sacrificing too much in the 
way of sensitivity to softer-spectrum and/or point- 
like sources, a set of parallel analyses was used. Each 
analysis is a variation on a common base analysis, 
optimized for better sensitivity to a particular type 
of source. In order to limit the number of addi- 
tional trials factors incurred, the number of paral- 
lel analyses was restricted to four. Two were opti- 
mized for sources with a Crab-like spectrum — one for 
point sources and the other for moderately extended 
(r = 0.2°) sources — and two more were optimized for 
hard-spectrum point-like and extended sources, using 
a reference source with a power-law spectrum and a 
spectral index of 2.0. 

The completely common elements of the analy- 
sis procedure consist of calibrating and cleaning the 
Cherenkov images and parameterizing them by second 
moments The technique used to stcreoscopically 
reconstruct the shower direction and impact parame- 
ter is likewise common to all analyses: however, im- 
ages used in stereoscopic reconstruction are required 
to exceed a minimum integrated charge {size) in dig- 
ital counts (dc), and the value of that requirement 
is analysis-specific as shown in Table HI Showers are 
reconstructed for all events in which at least three 
telescope images passed quality selection. 
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Table I Cuts Specific to Parallel Analyses 





Point Source 


Extended 
Source 


Standard 
Source 


size > 600dc 
(- 90p.e.) 
< 0.013 


size > 600dc 
< 0.055 


Hard Source 


size > lOOOdc 
(~ 150p.e.) 
< 0.013 


size > lOOOdc 
< 0.055 



Cosmic ray background is rejected using two means. 
The first is a pair of variables (mean scaled length 
(MSL) and mean scaled width (MSW) that summarize 
differences in image shape between gamma ray events 
and the majority of the cosmic ray background) @ and 
are apphed prior to generation of photon sky maps. 
The second is a cut on the square of the angular dis- 
tance (9^) between a reconstructed shower and the sky 
position of a potential source that is applied as part 
of that process. The cuts on mean scaled length and 
mean scaled width are common to all parallel analyses 
(0.05 < MSW < 1.06, 0.05 < MSL < 1.24), while the 
cut in 9^ is analysis-specific as shown in Table IH The 
residual cosmic ray background is estimated using the 
"ring-background" model 



4. Survey Detections 

Figure [2] shows the pre-trials significance at all 
points in the VERITAS Cygnus Region survey above 
76° galactic longitude, based on the standard ex- 
tended source analysis and including all data taken 
through Fall 2009. This region of the survey con- 
tains two clear detections of VHE gamma-ray emis- 
sion: the known source TeV J2032-|-4130[i^, and 
a newly discovered source, VER J2019-I-407, in the 
vicinity of the 7 Cygni SNR. As shown in Figure [H 
both sources are in proximity to a Fermi pulsar, and 
TeV J2032-h4130 is also coincident with the multi-TcV 
gamma-ray source MGRO 2031-1-41. 

In the full survey dataset, VER J2019+407 is de- 
tected at lOcr pre-trials and TeV J203-F4130 is de- 
tected at 5(T at the nominal position. Neither source 
exceeded 5cr pre-trials in the base survey. 



4.1. VER J201 9+407 

The base survey had shown a hint of extended emis- 
sion (at the 4(7 level) in the northwest corner of the 7 
Cygni SNR; survey-style follow-up data in Fall 2008 
and wobble observations taken in Spring 2009 raised 
the signal to the level of ~ 5 — 7cr pre-trials in the vari- 
ous survey analyses. In Fall 2009, a set of independent 
wobble observations were taken at 0.6° offsets about 
the position (RA,Dec) 20''19'"48", Dec = 40°54'00". 



The resulting significance map is shown in Figure 3(a)[ 
extended emission was detected with a maximum sig- 
nificance of 8.5 standard deviation {7.5a post-trials) 
within a search region of radius 0.25° around the 
central wobble position, using the standard extended 
source analysis from the survey. The analysis and 
search region were chosen a priori based on the the 
previous survey results. 

The centroid and intrinsic extension of the source 
are based on the Fall 2009 observations only, and 
were characterized by fitting an asymmetric two- 
dimensional Gaussian convolved with the point-spread 
function (PSF) of the instrument, to an acceptance- 
corrected uncorrclatcd excess map with 0.05° bins. 
The PSF was derived from a fit to Crab data taken 
in Fall 2009 and has a 68% containment radius of 
0.1°. We find a preliminary centroid position of RA = 
20'^19'"52.80'' ± 0.017°, Dec = 40°47'24.0" ± 0.023° 
and a preliminary extension of 0.16° ± 0.028° and 
0.11° ± 0.027° along the major and minor axes. 

The nature of VER J2019-f407 is stih unclear. As 
the centroid of VER J2019-I-407 is displaced from the 
Fermi pulsar position by over 0.5° (15pc for an as- 
sumed distance of 1.7kpcpT|). an association between 
the two sources seems unlikely. While it is not uncom- 
mon for VHE gamma-ray pulsar wind nebulae to be 
displaced from both the associated pulsar and the X- 
ray (or GcV gamma-ray) counterpart of the PWN, a 
displacement of this size would be rare and not easily 
justified. It seems more likely that the TcV emission is 
in some way associated with the northwestern portion 
of the SNR itself. It would be attractive to interpret 
the emission as being hadronic in origin and associ- 
ated with a shock-cloud interaction, but this would 
require sufficient density in HI to compensate for the 
relative lack of molecular material, particularly CO, 
at the VER J2019-I-407 position. There are two other 
puzzling features that must also be explained: the rea- 
son for the correlation between the TeV emission and 
the population of GeV electrons indicated by high- 
intensity regions in the radio contours at 1420 MHz 
and 4.85 GHz, and the confinement of emission to the 
northwestern portion of the remnant. One possible 
scenario is that the SNR was previously expanding 
within a low density bubble blown by winds from the 
progenitor star. In this case, the SNR has only re- 
cently begun to interact with relatively dense mate- 
rial in the interstellar medium, with the northwestern 
portion of the remnant being the first to encounter 
denser material. 

Uchiyama et al. have also argued for the presence 
of shock-cloud interactions in the northwestern part 
of the remnant on the basis of clumpy structures in 
0.7-lkcV ASCA satellite data, some of which are co- 
locatcd with VER J2019-h407[l3. However, the cor- 
responding identification of hard X-ray sources in the 
same region in the 4-7 keV band with shocked, high- 
density HI cloudlets rested on the identification of this 



eConf C091122 



4 



2009 Fermi Symposium, Washington, D.C., Nov. 2-5 




Galactic Longitude [Deg] 

Figure 2: Significance map of the portion of the VERITAS Cygnus Region Sky above 76° galactic longitude, including 
all follow-up data taken through Fall 2009. The error circles of two sources from the Fermi six- month bright source 
list[3], OFGL J2021. 5-1-4026 and OFGL J2032. 2-1-4122, both pulsars, are indicated by black circles, and the position 
and extent of MGRO J2031-|-4l|i] is indicated in blue. The position of TeV J2032-f4130 as measured by HEGRA and 
MAGIC is indicated by the dark and light blue crosses, respectivelvfisl] . 



region with EGRET source 3EG J2020+4017, which 
is now identified with OFGL J2021. 5+4026 and there- 
fore now localized well away from VER J2019+407. A 
Fermi upper limit on GeV gamma-ray emission within 
the VER J20194-407 emission region would therefore 
be essential to placing an upper limit on possible HI 
cloudlet density in the region. 



5. Base Survey Sensitivity and Limits 

A set of detailed survey simulations, coupled with 
observations of the Crab Nebula at multiple offsets, 
allows us to determine the a priori sensitivity of the 
survey analyses, not only to a point source with a 
Crab-like spectrum, but to harder-spectrum and/or 
significantly extended sources. In order to best re- 
produce the expected background conditions, blank 
(i.e. cosmic-ray dominated) survey fields were used to 
provide the background for these simulations. Blank 
fields at an appropriate range of zenith angles were 
pulled from the survey and arranged in a mocked-up 
'cell' of the survey grid around a test point. Show- 
ers that reconstruct at a distance greater than 1.7° 
from the center of the field of view are not used in 
the analysis of survey data. Therefore, only pointings 
where the center of the field of view is less than 1.7° 
from the test position are included in the simulation, 
as only these pointings contribute significantly to the 
sensitivity at the test point. 

Simulated gamma rays were then injected into each 



field in the simulated survey grid at the test position 
at a rate appropriate to the source spectrum and flux 
strength in question. In each case the gamma rays 
were simulated at an appropriate camera offset and 
matched as closely as possible to the background field 
in terms of zenith angle and azimuth. In the case of an 
extended source, the injection positions were smeared 
by a two-dimensional Gaussian with appropriate radii. 
The simulated grid is then analyzed using standard 
survey analysis procedures. 

Table II Survey Sensitivity 



Analysis 


So 

Spectral 
Index 


urce Properl 
Extension 
(radius) 


ies 
Flux 


Standard 
Point 
Standard 
Extended 


2.5(2.0) 
2.5(2.0) 


none 
0.2° 


4% Crab > 
200 GeV 
10% Crab 
> 200 GeV 


Hard Point 

Hard 
Extended 


2.0 
2.0 


none 
0.2° 


6.3% Crab 

> 500 GeV 
16% Crab 

> 500 GeV 



Preliminary results of the sensitivity studies, using a 
simulated survey grid of four-telescope array observa- 
tions taken at 20° zenith angle, have been completed. 
Using the results of theses studies, and accounting for 
variations in sensitivity due to array configuration and 
zenith angle, we derive preliminary values for the sen- 
sitivity of the various survey analyses to sources of 
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Figure 3: (a) Significance map of the region around VER J2019+407, based on Fall 2009 observations alone. 
Acceptance-corrected excess map of the region around VER J2019+407, using all data taken through Fall 2009. The 
black contours show the 1420MHz radio emission in the region (taken from the Canadian Galactic Plane Survey 
and clearly show the extent of the 7 Cygni remnant. The white-dashed line shows the lo" ellipse from the VER 
J2019+407 source extension fit; the circle is the error circle of OFGL J2021.5-f 4026 and the star is the position of the 
associated Fermi pulsar. 



particular spectral indices and extensions, as summa- 
rized in Table HTl Extrapolating from these sensitivi- 
ties and assuming a 3o" fluctuation, we calculate pre- 
liminary flux limits (for all points below 3(t in the 
survey region) of 3% of the Crab above 200 GeV for 
point sources and 8.5% of the Crab above 200 GeV 
for extended sources. 

The base Cygnus region survey sees no points at or 
above 5cr pre-trials. If we extrapolate from our sensi- 
tivity estimates and the population seen by the orig- 
inal H.E.S.S. surveyQ of the region between —30 < 
Z < 30, and correct for the area covered by the two sur- 
veys at or near full sensitivity, we would have expected 
to see between two and three sources in the base sur- 
vey. While these preliminary results are not conclu- 
sive, they do seem to indicate a difference in VHE 
gamma-ray source population between the Cygnus 
region and the Galactic center region surveyed by 
H.E.S.S. 



of the Crab Nebula for all points in the survey below 
3(7. Given our estimated sensitivity and the known 
source population from the original 2005 H.E.S.S. sur- 
vey of the Galactic Center region, we would have ex- 
pected to see 2-3 sources in the Cygnus region in the 
base survey at or above ha pre-trials, but see none — 
a result indicative of a population difference between 
the two regions. Follow-up observations within the 
survey region have also resulted in two source detec- 
tions: TeV J2032-I-4130 (which is hkely related to the 
coincident Fermi source OFGL J2032. 2+4122) and a 
newly-discovered region of extended emission, VER 
J019-f 407, overlapping the northwest portion of the 7 
Cygni SNR, which may or may not have a relationship 
to the nearby Fermi source OFGL J2021. 54-4026. 

VERITAS observations continue on areas of inter- 
est within the Cygnus Region survey. Studies of cor- 
relations between VERITAS and Fermi maps in this 
region are also planned. 



6. Summary and Outlook 

VERITAS has completed a 140-hour survey of a 
5 by 15 degree portion of the Cygnus region of the 
Galactic Plane, and simulation studies indicate that 
the base survey's average VHE point-source flux sensi- 
tivity is better than 4% of the Crab Nebula at energies 
above 200 GeV, and we derive preliminary flux limits 
of 3% of the Crab Nebula for point sources and 8.5% 
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